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Abstract 

We propose an discontinuous Galerkin local orthogonal decomposi¬ 
tion multiscale method for convection-diffusion problems with rough, 
heterogeneous, and highly varying coefficients. The properties of the 
multiscale method and the discontinuous Galerkin method allows us 
to better cope with multiscale features as well as interior/boundary 
layers in the solution. In the proposed method the trail and test spaces 
are spanned by a corrected basis computed on localized patches of size 
0{H log{H~^)), where H is the mesh size. We prove convergence rates 
independent of the variation in the coefficients and present numerical 
experiments which verify the analytical findings. 


1 Introduction 

In this paper we consider numerical approximation of convection-diffusion 
problems with possible strong convection and with rough, heterogeneous, and 
highly varying coefficients, without assumption on scale separation or period¬ 
icity. This class of problems, normally refereed to as multiscale problem, are 
know to be very computational demanding and arise in many different areas 
of the engineering sciences, e.g., in porous media ffow and composite materi¬ 
als. More precisely, we consider the following convection-diffusion equation; 
given any / G we seek u G = {n G H^{Q) \ n|r = 0} such that 

— V ■ AVu + b ■ Vu = f in G, (1) 

is fulhlled in a weak sense, where G is the computational domain with bound¬ 
ary r. The multiscale coefficients A, b will be specihed later. There are two 
key issues which make classical conforming hnite element methods perform 
badly for these kind of problems, 

• the multiscale features of the coefficient need to be resolved by the 
hnite element mesh and 
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• strong convection leads to boundary and interior layers in the solution 
which need to be resolved. 

To overcome the lack of performance using classical hnite element meth¬ 
ods in the case of multiscale features in the coefficient many different so called 
multiscale methods have been proposed, see [251 ESI ESI 0 ESI 1121 EDI ET] 
among others, which perform localized hne scale computations to construct 
a different basis or a modihed coarse scale operator. Common to the afore¬ 
mentioned approaches is that the performance of the method rely strongly 
on scale separation or periodicity of the diffusion coefficients. There is also 
approaches which perform well without scale separation or periodicity in the 
diffusion coefficient but to high computational cost by either having to solve 
eigenvalue problems [2] or where the support of the localized patches is large 
[371 H] . See also [38] . 

In the variational multiscale method (VMS) framework [251 ED] the so¬ 
lution space is split into coarse and hne scale contribution. This idea was 
employed for multiscale problems in a adaptive setting for classical hnite 
element in [an [311 [32] and to the discontinuous Galerkin (DG) method in 
ig. A further development is the local orthogonal decomposition (LOD) 
method, see [351 EDI ED] for classical hnite element and [15] for DG meth¬ 
ods. The LOD operates in linear complexity without any assumptions on 
scale separation or periodicity and the trail and test spaces are spanned by 
a corrected basis function computed on patches of size 0{H\og{H~^)). The 
LOD has e.g. been applied to eigenvalue problems [35], non-linear elliptic 
problems [21], non-linear Schrodinger equation m, and in Petrov-Galerkin 
formulation [TB] . 

There is a vast literature on numerical methods for convection dominated 
problems, we reefer to [281E31E7], among others. There has also been a lot 
of work on DG methods, we refer to [391 [331 [31E9] for some early work and to 
[HI [221 iQl [9] and references therein for recent development and a literature 
review. DG methods exhibit attractive properties for convection dominated 
problems, e.g., they have enhanced stability properties, good conservation 
property of the state variable, and the use of complex and/or irregular meshes 
are admissible. For multiscale methods for convection-dihusion problems, see 
e.g. mSDEH]. 

In this paper we extended the analysis of the discontinuous Galerkin local 
orthogonal decomposition (DG-LOD) [T5] to convection-diffusion problems. 
For problems with strong convection using the standard LOD won’t suffice, 
since convergence can no longer be guarantied. Instead we propose to include 
the convective term in the computations of the corrected basis functions. We 
prove convergence results under some assumptions of the magnitude of the 
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convection and present a series of numerical experiment to verify the analytic 
findings. For problems with weak convection it is not necessary to include 
the convective part 0. 

The outline of this paper is as follows. In section |2] the discrete setting 
and underlying DG method is presented. In section [3] the multiscale decom¬ 
position, the DG-LOD, and the corresponding convergence result are stated. 
In Section 0] numerical experiments are presented. Finally, the proofs for 
some of the theoretical results are given in Section [5l 


2 Preliminaries 

In this section we present some notations and properties frequently used in 
the paper. 


2.1 Setting 

Let C for d = 2, 3 be a polygonal domain with Lipschitz boundary F. 
We assume that: the diffusion coefficients, A G has uniform 

spectral bounds 0 < a,/? < oo, dehned by 


0 < a := ess inf inf 


(A(x)v) ■ V 


xen DeM‘^\{o} v ■ v 


< ess sup sup 


(y4(x)n) 




=: (3 <oo, (2) 


and the convective coefficient, b G is divergence free 


V ■ b(a:) = 0 a.e. x G 12. (3) 

We denote Ca = (/9/a)^/^. 

We will consider a coarse and a hne mesh, with mesh function h and H 
respectively. Furthermore, we assume that the hne mesh resolve and that 
the coarse mesh do not resolve the hne scale features in the coefficients. Let 
Tfc, for k = {h, id}, denote a shape-regular subdivision of 12 into (closed) 
regular simplexes or into quadrilaterals/hexahedras (d = 2/d = 3), given a 
mesh function fc : 7/ —)■ M dehned as k := diam(T) G Po{7k) for all T G Tk- 
Also, let VfcU denote the 7fc-broken gradient dehned as (Vn)|r = VuIt for all 
T EPk- For simplicity we will also assume that Tk is conforming in the sense 
that no hanging nodes are allowed, but the analysis can easily be extend to 
non-conforming meshes with a hnite number of hanging nodes on each edge. 
Let T be the reference simplex or (hyper)cube. We dehne Vp{T) to be the 
space of polynomials of degree less than or equal to p if T is a simplex, or the 
space of polynomials of degree less than or equal to p, in each variable, if T 
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is a (hyper)cube. The space of discontinuous piecewise polynomial function 
is defined by 

Pp{Tk) := {n : ^ M I VT e Tfc, v\t o Ft G Pp(f)}, (4) 

where TV : T ^ T, T G 7^ is a family of element maps. We will work 
with the spaces Vk '■= Pi(Tk). Let UpiTk) ■ LF‘{Vt) —)■ Pp{Tk) denote the 
projection onto PpiTk)- Also, let £k denote the set of all edges in Tk where 
£k{£l) and Tfc(r) denote the set of interior and boundary edges, respectively. 
Given that and T~ are two adjacent elements in 7^ sharing an edge 
e = n T~ G £k{£t), let Vf. be the the unit normal vector pointing from 
T~ to T"*", and for e G Tfc(r) let Ve be outward unit normal of G. For any 
V G PpiTk) we denote the value on edge e G ^(G) as = n|enT±- The jump 
and average of n G PpiTk) is dehned as, [n] = v~ —v^ and {n} = (n“ + v^)/2 
respectively for e G and [n] = {n} = v\e for e G Tfc(r). For a real 

number x we dehne its negative part as x® = l/2(|a;| — x). 

Let 0 < C* < cxD denote any generic constant that neither depends on the 
mesh size or the variables A and b; then a < b abbreviates the inequality 
a < Cb. 

2.2 Discontinuous Galerkin discretization 

For simplicity let the bilinear form ah{-,-) : x Vh —?■ M, given any mesh 

function h : Q ^ Po(7h), be split into two parts 

ah{u,v):=a‘l^{u,v) + al{u,v), ( 5 ) 

where a^(-, •) represents the diffusion part and a/j(-, •) represents the convec¬ 
tion part. The diffusion part is approximated using a symmetric interior 
penalty method 

atiu,v) := {AThU,Thv)LHn) + 

- {{Ue ■ ATu}, [v])L^e) “ ' AVv}, ML2(e))) , 

where cXe is a constant, depending on the diffusion, large enough to make 
ajj(-, •) coercive. The convective part is approximated by 

ah{u,v) := ih-ThU,v)L2(n)+ ^ ibe[u],[v])L2(e) 

- • b[M],{n})i 2 (e) + ((z/e • b)®M,'i;)i 2 (e), 

eG5^(r) 
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where upwind is imposed choosing the stabilization term as 6e = |b ■ i^e|/2 

0 . 

The following dehnitions and results are needed both on the hne and 
coarse scale, for this sake let k = {h, if}. The energy norm on Vk is given by 


e€£k 

eGffc 


lL2(e), 


( 8 ) 


I 1112 111 1112 I 111 1112 

Tlllfc = IItIIIm + llTlIlfc.c- 


From Theorem 2.2 in [30] we have that for each v G Vk, there exist an 
averaging operator —?■ Vfc fl H^{kl) with the following property 


|Vt(« - + II v‘(« - < 5^ TllHlli.,.,. (9) 


eGSk 


In the error analysis we will also need a localized energy norm, dehned in a 
domain a; C fi (aligned with the mesh Tk) as 


I /c,d,cj 




e€£k 

enco^O 


\ k,c,u} 


E PI'MWIh.), 


\ k.LJ 


e&£k 

eno^O 

_ I I l-y I ' 


I k.d.LJ 


+ k' 


Ifc.C.CJ* 


( 10 ) 


3 Multiscale method 


In this section we preset the multiscale decomposition and extend the results 
in na to convection-diffusion problems. For the constants in the convergence 
results to be stable we assume the following relation of the convective term 


^/liybji^' 




a 


( 11 ) 


How the magnitude of flTT]) affects the convergence of the method is investi¬ 
gated in the numerical experiments. 
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3.1 Multiscale decomposition 

In order to do the multiscale decomposition the problem is divided into a 
coarse and a hne scale. To this end let Th and 7^, with the respective 
mesh function H and h, denote the two different subdivisions, where 77 is 
constructed by some (possible adaptive) rehnements of Th- 

The aim of this section is to construct a coarse hnite element space based 
on Th, which takes the hne scale behavior of the data into account. We 
assume that the mesh 77 resolves the variation in the data, i.e., the solution 
to: hnd Uh E Vh such that 

ah(uh,v) = F(v) for all e V/i, (12) 

gives a sufficiently good approximation of the weak solution u to o. Note 
however that never have to be computed in practice, it only acts as a 
reference solution. We introduce a coarse projection operator Ilji/ ;= ni(77/) 
and let the hne scale reminder space be dehned by the kernel of Ilji/, i.e., 

yf := {vEVh\ HhV = 0} C Vh- (13) 

The coarse projection operator has the following approximation and stability 
properties. 

Lemma 1. For any v eVh and T G Th, the approximation property 

H\t^\\v - UhvWl^t) < a~^^'^\\\v\\\h,T, (14) 

and stability estimate 

ll|nH^|||rr<C.|||n|||h, (15) 

is satisfied, with 

( 16 ) 

Proof. The approximation property follows directly from m Lemma 5]. Let 
Ch ■ ^ n Vh be a Clement type interpolation operator proposed in 

[HI Section 6] which satisfy 

\\TChu\\l2{t) + - Chu)\\l'2{t) < II V-u||2,2(^^), (17) 

where uf = int(U{T' ETh I T n T 7 ^ 0}) are the union of all elements that 
share a edge with T. We dehne the conforming function Vc = CnTfiv using 
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averaging operator in ([2]). We obtain 

TgTh 

+ E (^IIK-nHIlWi + llst^l-c-nHir) , 18 ) 

eGShii^'JT'o) 

~ ^ ~ ^ovWh^T) + (jp + llwn - n||i2(r)^ 

using that Ho ;= no(TH-) is the L^-projection onto constants, a trace inequal¬ 
ity, and stability of IIj:/. Next, using that 

WChT^v - v\\L2^n) < \\Ch^v -X^^v\\L2^n) + Wl^v - v\\L2(n) 

< H\\VT^v\\L2(n) + \\XhV - v\\L2(n) (19) 

< a-X^H\\\v\\\i, 

in ffT8|l concludes the proof. □ 

The following lemma shows that for every vh G Vh there exist a (non¬ 
unique) V G HJj^Vh G Vh in the preimage of If// which is conforming. 

Lemma 2. For each vh G Vh, there exist a v E Vh V\ H^{Q) such that 
Uhv = Vh, lll^^lllh < Ca\\\vh\\\h, and supp(n) C supp(X^nH)- 

Proof. Follows directly from [151 Lemma 6], since v G H^{Q). □ 

The next step is to split any v E Vh into some coarse part based on Th, 
such that the hne scale reminder in the space is sufhciently small. A naive 
way to do this splitting is to use a L^-orthogonal split. An alternative dehni- 
tion of the coarse space is Vh = ^nVh- A set of basis functions that span Vh 
is the element-wise Lagrange basis functions {Atj | T E Th, j = 1, • • • 
where r = {1 + d) for simplexes or r = 2'^ for quadrilaterals/hexahedra. The 
space Vh is known to give poor approximation properties if Th does not re¬ 
solve the variable coefficients in ([1]). We will use another choice, see [361 [T5] . 
based on a/i(-, •), to construct a space of corrected basis functions. To this 
end, we dehne a hne scale projection operator : V/i —?■ by 

ah{^v,w) = ah{v,w) for all w G (20) 

and let the corrected coarse space be dehned as 

Vf}^ ■.= {!-d)VH. ( 21 ) 
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The corrected space are spanned by corrected basis functions := {Xtj — 
<pT,j I T G Th, j = ^, ■ ■ ■ which can be computed as: for all T G Th, j = 
1,... ,r hnd (pT,] ^ such that 

ah{4>T,j, v) = ah{XT,j, v) for all v G (22) 

Note that, dim(V)^^) = dim(VH)- From fl2l|) we have that any Vh G Vh can 
be decomposed into a coarse G and a hne G scale contribution, 
Vh = v]}^ + vL 

Lemma 3 (Stability of the corrected basis function). For all T G Th, j = 
1 , ... ,r, the following estimate 

|||0T,h - Arjlllh < C^lT^‘^\\H~^XT,j\\L-^{n) (23) 

holds, where C(j, = {C\ + 

Proof. Let v = Xtj — bxj G where brj G Hq(T), Unbrj = Xtj, 
III^Tjlllh < IIAtjII|h from Lemma[2l We have 

I \\4>T,h — AtjI I Ih ^ Clh{4>T,h — Xtj, 4>T,h — Xtj) 

= <Vh{4>T,h — Xtj, V — Xtj) = cihi4>T,h — Xtj, bTj) (24) 

= ah{4'T,h - Xtj, bTj) + (b ■ Vh{(j)T,h - Xtj), bTj)Ti{n) ■ 

Using that the diffusion part in fl2T)) of the bilinear form is continuous in 
(Vh X Vh) with the constant Ca, Lemma [2l and a inverse inequality, we get 

(4i{4>T,h — Xtj, bTj) < C'yl|||0T,h — -^TjIlIhlll&Tjlllh 

~ CW\\(t>T,h - Arj|||h|||Arj|||H (25) 

< C\j3^/‘^\\\(i)T,h — -^Tjlllhll-f^ ^Xtj\\l'^{T)- 

For the convection part in fl241) . we have 
(b ■ Th{4>T,h — Xtj), &Tj)l2(o) 

< \\Hh ■'Vh{4>T,h — XTj)\\L‘2{n)\\H (26) 

< ||hfb||i,oo(Q)||Vh(0T,h — ATj)||L2(o)||Ff ^XTj\\L^(n), 

and obtain 

\\\4>T,h — ATjlllh < ^Atj||l2(q). 

with Cj = {C\ + ||F/'b||Lcx>(Q)Q:"i)F2, 


(27) 

□ 


3.2 Ideal discontinuous Galerkin multiscale method 

An ideal multiscale method seeks G such that 

ah{u^^, v) = F{v) for all v e (28) 

Note that, to construct in the space a variational problem has to be 
solved on the whole domain hi for each basis function, which is not feasible 
for real computations. The following theorem shows the convergence of the 
ideal (non-realistic) multiscale method. 

Theorem 4. Let Uh G Vh he the solution to flT^ . and G be the 
solution to (l28|) . then 

11!«/. - “SI 11 ;S - n„/)| ( 29 ) 

holds, with Cl = Ca + ||hfb|hoo(Q)a“^ 

Proof. See Section O □ 

3.3 Discontinuous Galerkin multiscale method 

The fast decay of the corrected basis functions (Lemma [6]), motivates us to 
solve the corrector functions on localized patches. This introduces a local¬ 
ization error, but choosing the patch size as 0{Hlog{H~^)) (Theorem [7]) 
the localization error has the same convergence rate as the ideal multiscale 
method in Theorem 01 The corrector functions are solved on element patches, 
defined as follows. 

Definition 5. For all T G Th, let be a patch centered around element T 
with size L, defined as 

uf := int(T), 

(30) 

cjy := int(U{T'G I T 0}), L = l,2,.... 

See Figure]^ for an illustration. 

The localized corrector functions are calculated as follows: for all {T G 
Tff, j = 1, • • •, r} hud G = {n G yf I = 0} such that 

a/j(0^j,n) = a/j(ATy, w), for all n G 12^(0;^). (31) 

The decay of the corrected basis function is given in the following lemma. 
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Figure 1: Example of a patch with one layer, two layers and three 
layers cafi, centered around element T. 

Lemma 6 . For all T G Th, j = 1, • • • where (pTj is the solution to fl2^ 
and is the solution to fl36|) . the following estimate 

I I I ^T,j - 4>T,j I I L ~ ^27^ 111 i^T,j - 4>T,j Ilk (32) 

£(fc-l)-l 

holds, where L = ik is the size of the patch, 0 < 7 = 2 ^fcn+i) < 

C 2 = CcC(^{l + CaCs), and C 3 = C(C^ + ||Frb||ioo(Q)a“^), where C is a 
generic constants neither depending on the mesh size, the size of the patches, 
or the problem data. 

Proof. See Section |5l □ 

The space of localized corrected basis function is dehned by ; = 

{<pT,j ~ ^T,j I T G Th, r = 1 ,..., r}. The DG multiscale method reads: hnd 

ms’,L yms,L , , , , 

t ^fj oU-Cll LllcLL 

ah{uff^'^,v) = F{v) for all n G (33) 

An error bound for the DG multiscale method using a localized corrected 
basis is given in Theorem [3 Note that it is only the hrst term |||m — M/i|||h 
in Theorem [3 that depends on the regularity of u. 

Theorem 7. Let Uh G Vh and G be the solutions to isa and 

(I31?]) . respectively. Then 

11 |w - w^’^l I Ih < 11 |w - wk I k + C,a-^/^\\H{f - Unf) 

+ G5||i7-iLoo(f,)L''/Vll/l|L2(0) 
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holds, where L is the size of the patches, Ci is a constant defined in Theo- 
rem\^ 0 < j < 1 andC^ = where C 4 = + CaCsY is 

defined in Lemma fTl and C 2 and 7 are defined in Lemma\B 

Proof. See Section O □ 

Remark 8. Theorem\^is simplified to, 

|||m - < |||m - M/i|||h + Ci\\H\\l^^q). (35) 

given that the patch size is chosen as L = [C'log(R“^)] with an appropriate 
C and 11 / 112,2 = 1. In the numerical experiments we choose C = 2. 

Remark 9. If the convective term is small it is not necessary to include it in 
the computation of the correctors JH]/. Instead the following correctors can 
be used: for all {T G Tff,/ = 1,..., r} find E such that 

v), for all v E V^{u^). (36) 

This gives the right convergence results if 

O = 1 (37) 

compared to ffTT]) if the convective term is included. 

4 Numerical experiment 

We consider the domain ^2 = [0,1] x [0,1] and the forcing function / = 
1 + cos(27ra:) cos(27r|/). The localization parameter which determine the size 
of the patches is chosen as L = [21og(R“^)], i.e., the size of the patches are 
2 Rlog(R“^). Consider a coarse quadrilateral mesh, Th, of size H = 2 ~\ i = 
2, 3,4, 5. The corrector functions are solved on sub-grids of the quadrilateral 
mesh, Th, where h = 2“^. We consider three different permeabilities: Ai = 1, 
A 2 = ^ 2 ( 2 /) which is piecewise constant with respect to a Cartesian grid of 
width 2 “® in y-direction taking the values 1 or 0 . 01 , and A^ = Afix, y) which 
is piecewise constant with respect to a Cartesian grid of width 2“® both in the 
X- and y-directions, bounded below by a = 0.05 and has a maximum ratio 
fi/a = 4 • 10^. The permeability A^ is taken from the 31 layer in the SPE 
10 benchmark problem, see http:www.spe.org/web/csp/. The diffusion 
coefficients A 2 and ^3 are illustrated in Figure Hi For the convection term 
we consider: b = [C, 0 ], for different values of C. 
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Figure 2: The diffusion coefficients A 2 and in log scale. 



Figure 3: The number degrees of freedom (iV^o/) vs. the relative error in 
energy-norm, for different sizes of the convection term, C. 


To investigate how the error in relative energy-norm, 11 \uh—u^^’^\\\/\ \ |Mh| 11, 
depends on the magnitude of the convection we consider: Ai and b = [C, 0] 
with C = {32, 64,128}. Figure [3] shows the convergence in energy-norm as a 
function of the coarse mesh size H for the different values of C. 

Also, to see the effect of heterogeneous diffusion of the error in the relative 
energy-norm, \\\uh — 11/|||m/i| 11, we consider: Figure 0] which shows the 

error in relative energy-norm using A 2 and b = [1,0] and Figure 0] which 
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shows the error in relative energy-norm using and b = [512, 0]. 



Figure 4: The number degrees of freedom {Ndof) vs. the relative error in 
energy-norm, using a high contrast diffusion coefficients A 2 and b = [1,0]. 
The dotted line corresponds to . 

We obtain convergence of the DG multiscale method to a reference 
solution in the relative energy-norm, \\\uh — 11 /|||m/i|11 , independent of 

the variation in the coefficients or regularity of the underlying solution. 


5 Proofs from Section [3] 

In this section we state the proofs of the main results which was postponed 
from in section [3l To this end we start by proving some technical lemmas in 
Section 15.11 which we use to prove the main results in Section 15.21 

5.1 Technical lemmas 

In the proofs of the main results. Theorem 01 Lemma [ 6 l and Theorem [71 we 
will need some definitions and technical lemmas stated below. 

Continuity of the DG bilinear form for convective problems can be proven 
on a orthogonal subset of Vh- The space is an orthogonal subset of Vh but 
on a coarse scale. 
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Figure 5: The number degrees of freedom {N^of) vs. the relative error in 
energy-norm, using a high contrast diffusion coefficients A 3 and b = [512, 0]. 
The dotted line corresponds to . 


Lemma 10 (Continuity in {Vh x V^) and (V^ x V^)). For all, {u, n) e x Vh 
or in V/i X , it holds 


<Cr 




(38) 


where 

a = C'^ + ||i/b|Ucc(o)«-h (39) 

Proof. Since is continuous in (Vh x Vh) with the constant Ca, continuity 
in (V^ X Vh) follows from V^ C Vh- For the convective part a)), we have 


a^(v,w) = ^ (b • Vv,w)l^T) + ^iK[v], M)L2(e) 

TGTIi 

,w)L2(e) 

eG£k{Q) eG5fc(r) 

~ (||b||L°°(r)||Vn||L2(r)||u'- nH'u;||L2(T)) 

T&Th 
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where S~^,S E Th and e = S~^ H S . Using a discrete Cauchy-Schwartz 
inequality and summing over the coarse elements, we get 

a''{v,w) < a~^^‘^\\Hh\\Loo^n)\\\v\mH~\w - 


< lli^bl 


L^(Q,)Ct 


- 1 | 


PlIIhllFlIlh, 


(41) 


which concludes the proof for (V^ x V^). The proof of (V^ x Vh) is obtained 
by hrst integrating (b ■ 'Vu,v)l^(^t) by parts. □ 

The following cut-off function will be frequently used in the proof of the 
main results. 

Definition 11. The function G Poilh), for D > d, is a cut off function 
fulfilling the following condition 


c 


d,Di _ 


T 


= 1 , 


_ PI 

ST — U, 


(42) 


\\[CT^]\\L-°{ShiT)) 


< 


L°°{T) 


{D-d)H\T^ 
and = 0, for allT e Th- 

For the cut off function has the following stability property. 

Lemma 12. For any v eVh and from Definition [771 the estimate, 

^d,D 


IllCr’ ^lllh < 

holds, where C(^ = {C\ + ||hb|| 2 ,oo(o)/a)^/^. 

Proof. For the diffusion part we use the following result from [15 

111(1-C?>llkd<c'^lll^l 

and focus on the convective part. We obtain 

111(1-d>iiiL 


(43) 


I h,Q\iJ, 


L-1 

T 


(44) 


eeSh 

< E (|I'>FI 


lL2(e) 


< 


=°(S+U5-)II-^ 


e^Sh- 


( 45 ) 


< 


T^Th-. 

't' 

||hb||j;,oo(Q) 


eDuj^ ^7^0 


a 
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using [vw] = {u}[w] + {w}[u], the triangle inequality, and a trace inequality. 
The proof is concluded using flH|) and (USD . □ 


The following lemmas will be necessary in order to prove Theorem [3 
Lemma 13. The following estimate, 

III ^i(0T,i-0T,j)lllh ^ Y l^irill<^Ti-</>rjlllL (46) 

TGTH,j=l,--;r T€TH,j=l, — ,r 

holds, where C 4 = C^Cf{l + CaCs)^- 

Proof. The proof is analogous with the proof of Lemma 12 in [15]. □ 

5.2 Proof of main results 

We are now ready to prove, Theorem 01 Lemma [HI and Theorem [3 

Theorem Let us decompose Uh into a coarse contribution, G and 
a hne scale remainder, G V^, i.e., Uh = . For we have 

llk-^lllh ~ ahiv^,v^) = ahiuh,v^) = if,v^)L2^^-) 

= (/ - Unf, - UHV^)L^(^n) 

< \\H{f - UHf)\\mn)\\H-\vf - nHvf)\\L2^n) 
<a-^/^H{f-IlHf)\\LHn)\\\v^\\\w 


(47) 


Using continuity, we get 

\\\uh- 11 Ih < Muh - UH - n^) = dhiuH - uh - v^) 

<a|||w/^-wS1llh|||wh-^nilh, 

which concludes the proof together with fiT7|) - 


(48) 


□ 


Lemma 0 Dehne e := (pTj — 4>T,j where cfxj G and <p!f j G {ujff). We 
have 


< 


ah(e, 0Tj - 0r 7 ) = «/i(e, (frj - v) < Ud | |e| | |h| | I^Tj - f 11 |h- (49) 


Furthermore from Lemma [21 there exist a n = Cr — hr G {ojf) such 

that n^^hr = T\-H{CY’^(t>T,j) and IHhrlHh < CA\\\P-H{CY'^(t^T,j)\\\H, we have 


kllL ^ C'c (^111(1 - Ct ^’'^)<^Tj|||h + |||hr|||h) 


(50) 
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where 


'^Ca\\\^hCt ^'^(t>T,j\\\H — Ca\\\^hO- — Ct ^'^)(t>T,j\\\H 
< c^c.iiKi - c^-'’^)0T,lllh < c'aC'AIII0t,IIL,^\4-i 


using Lemma [2|, Lemma [T], and Lemma [T21 We obtain, 


(51) 


< 


C 2 ||| 0 ^ 


Tjl Nh,ov^“i’ 


(52) 


where = CcC((l + CaCs) from fl50l) and flHTl) . 

The next step in the proof is to construct a recursive relation which will 
be used to prove the decay of the correctors. To this end, let ik = L — 1, 
and dehne another the cut off function, rjjp := — and 

the patch ujtp := for m = 0,1,..., \lk/ (£ + 1) — Ij. Note that 

C cj™. We obtain 


111 4>T,j 111 h,Q\^jp <111 hT 111 h < Oft {r]^c()T,j , -^T 0T,i) • 


(53) 


To shorten the proof we refer to the following inequality 

a'^(?7™0T,i,'^T0Tj) < a'^(0Tj, (hT)^0Tj - br) + —^\\\4>t, 


'' 'h,u)'!p\u), 




. (54) 


where (?7 t )^0T,i — £ V^, in the proof of Lemma 10 in [15]. We focus on 

the convection term, since the cut of function is piecewise constant it follows 
that 

(b • Vr7™0Tj, VT^T,j)L2{s) = (b ■ (hr )^<(’Tj)l2 (5) (55) 

for all S E Th- Using the following equalities from (Appendix A in |15j ) 

{vw}[vw] = {w}[v^w\ — Hityjlnjltt;} + l/4[n]{n}[t(;] [to], 

[vw][vw] = [w][v^w] - 1/A[v]^[w]^ + [v]^wy, 


and fl5^ . we obtain 

a^(?7™0Tj,hr0rj) = a''(0rj, (hT)^0Tj) 

+ {i^e-h[r]tp]{(j)T,j},{VT}{(pT,j})me) 

ee£hi^) (57) 

- l/4(i/e • h[ritp]{(PTj}, {hr }[0T,i])L2(e) 

- l/ 4 (&e[h^]^ [ 0 T,jf )L 2 (e) + ( 6 e[hT]^ {(J^TjV)L^e))- 
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The sum over the edges terms can be bounded using that || [^r] IU°°(t) ^ 
\\h\\L^(T)/ H\t, \\{VT}\\L-°{n) < 1, \\h\\L^^T)/H\Ti < 1, and a trace inequality. 
We obtain 


eGf/i(r 2 ): 

+ II II 


< 




H b||ioo(e) 


1111 1,2(7^+ 

UT-) 


en(w 5 ?\i:D 5 ?)^o 


< 


E 


\Hh\ 


L-^iT) II „_i 


H (0TJ - ^H4’T,j)\\L2{T) 


T&Th- 

Tn{urp\ijrp)^o 

1112 


£a 


Td\\\h,iu^\si:p+^)- 


(58) 


Combining the results, we have 


lll0Tj|||ft,o\s5- < a(0Tj, iVTY<pT,j - hr) + a{(j)T,j,bT) 


+ i 


-1 


Cl + 


Ilifbll 


L°“( 0 ) 


a 


III'#'' 


(59) 


'TilNh,(C5j‘\i:5^+b’ 


where br has support in \ such that {j}'!^Y(t>T,j — and 

III^TlIlh < CA\\\^H{{'nTf(t)T,j)\\\H, see LemmaH We have 

— bx) = 0. (60) 

For all T G Tb the operator IIj:/ is stable in the L^(T)-norm, we have 






= Cl {mHiiv^nTMiH + mHiivmTMiiH) ■ 

For the hrst term in flOTH we refer to the result 

^m\2 \ll |2 ^A\ 


(61) 


|nH((t)”) ^j)\\\d,H S: ■^lpTj|||j^2„,j,y+i, 


(62) 
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from [15] and for the second them we have 


ee£H{n) 

= ||hf ^^\\l°°{T)\\{Vt)‘^ ~ ^oiVT)‘^\\‘L<=°{T)\\^T,j — 
TgTh 

^ ||hfb||Loc(r) 2 

~ q,£2 


( 63 ) 


We obtain 


\j\\\h,n\ujjp- 


= C3r^\\\(j)T,j\\\h,n\uii^+^ 


^ p -1 I ^2 I ll^b||ioo(r)\ 2 


a 


(64) 


where C 3 = C{C\ + ||-ffb||ioo(Q)CK and C is the generic constant hidden in 
We have 

for any m = 0 , 1 ,..., \fk/{i + 1 )J — 1 , which we can use recursively as 


Ill 0 r,, 


J 111 h,Q\L0j 


< 




J 111 h,Q,\u}^ 




|2 


( 66 ) 


Note that /c/2 is a lower bound of ik/{E + 1). Equation together with 
dSSj), gives 


lll'/’Tj - 0hlllh ^ _ Arjillh- 

which concludes the proof is concluded. 

Theorem^ Using the triangle inequality, we have 


\u — u 


/ms,L I 
H \ 


< \ \\u - Uh\\\h + \ \\Uh - U 


ms,L I 
H \ 


h- 


(67) 

□ 

( 68 ) 


Note that, Uh G V/i, can be decomposed into a coarse, G and a hne, 
G scale contribution, i.e., Uh = ■ Also, let G be 

chosen such that Uhv^^’^ = Uhv^- We have 

\\\Uh-Ujj’ \\\h^ah[Uh-u^ ’ ,Uh-UH ) 

K ms.L ms,L\ /'m\ 

-Ufj’ ,Uh-Vjj ’ ) (69) 


< 


Cc\\\Uh - U^''’^\\\h\\\Uh - V 


ms,L I 
H \ 
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and obtain 


111 ms.L III ^ 111 111 

\\\u-Ujj' lllh <|||m - Mhlllh 

+ C'c(|||wh-^rillh + |||n^^-nr"'ll|h). 

The first term in flTOj) implies that the reference mesh need to be sufficiently 
hne to get a sufficient approximation. The second term is approximated 
using flT7|) . i.e. 


IIK - t>rillh < - n„)/||i,(n,, 

and for the last term in we have, 


(71) 


I,,ms _ 7 , 


ms,L 111 2 _ 

Ih ~ 


'^H^Ti^j){^T,h 0T,j)lllh 

TeTH,j=i, — ,r 

<C,L^ \vH:TixM\^T,H - <Ph\\\l 

T&TH,j=^,--;r 


T&TH,j=^,-,r 

using Lemma dSl and Lemma O 
We obtain, using Lemma [3l that 

l^ir,r(2^i)nil0T,h - Arjilih 

TeTH,j=l,-,r 

<cl Y1 \\H-'vs‘T(^t)\Tj\\hm 

TeTH,j=i,—,r 


(72) 


< ClPW Y, H^''>H:T(x,)>.Tj\\h(n, 

T(iTH,j=l,-,r ( 73 ) 

^Cln Y ff“‘»«.Tfe)nH(AT,,-.#.Tj)|li.(n| 

T’GTTt, j=l,.--, 7 ’ 

holds and we conclude the proof. 

□ 
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